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A  computer-simulation  model  is  used  to  study  the  percolation  phase  transition  ot  stiff 
ehains/sticks  on  a  square  lattice.  Ry  varying  the  chain  length,  effects  of  the  chain  length  {L  ) 
on  the  percolation  transition  is  explored.  Ihe  percolation  exponents  seem  to  remain 
unchanged  with  L  .  The  percolation  threshold  />,  depends  strongly  tm  the  chain-length,  and 
shows  a  power-law  dependence  p  {I,  )  ~  I.  '  '. 


1.  Introduction 


Because  of  its  diverse  applicabilities  particularly  in  fluid  flow  through  porous 
media  (i.c.  marine  sediments),  percolation  phenomena  have  been  the  subject 
of  active  research  for  a  long  time  |1.2|.  Site  and  bond  percolations  arc  mostly 
studied,  and  have  been  quite  suecessful  in  understanding  the  physical  and 
chemical  properties  of  inhomogeneous  materials  (alloys  in  particular).  In  most 
of  the  polymeric  materials  i.e.  the  gel  networks,  the  basic  units/segments  that 
form  the  incipient  infinite  network,  are  chains  and  loops  of  various  sizes  [3].  In 
order  to  understand  the  strength  of  connectivity  at  the  onset  of  the  infinite 
network,  it  would  be  interesting  to  study  a  model  of  percolating  networks  with 
chains  and  loops  of  various  sizes  as  basic  units.  Site  and  bond  percolation 
mechanisms  may  not  be  adequate  to  describe  the  conformational  properties  ot 
such  random  networks.  On  the  other  hand,  it  is  rather  difficult  to  consider  both 
loops  and  chains  of  various  sizes  and  distributions,  because  of  limitations  on 
computer  memory  and  time.  Therefore,  it  might  be  useful  to  investigate  the 
percolation  of  chains. 

Drory  et  al.  [4|  have  recently  studied  the  percolation  threshold  of  permeable 
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objects  ot  various  shapes  in  continuum  space.  They  tound  that  the  percolation 
threshold  depends  strongly  on  the  size  and  shapes  (i.e.  aspect  ratii)  tor 
elongated  boxes,  radius  tor  spheres,  etc).  Here  we  have  investigated  the 
percolation  ot  slitt-chains  (sticks  or  thin  rods)  ot  various  sizes  in  a  two- 
dimensional  discrete  lattice:  we  also  find  a  strong  dependence  ot  the  percola¬ 
tion  threshold  on  the  chain  length.  Such  studies  may  be  helpful  in  short  tibres 
epoxy  resin  composites  |5|.  A  related  study  on  jamming  coverage  has  been 
recently  carried  out  by  Svrakic  and  Henkel  (b|  for  an  irreversible  deposition  of 
mixtures  of  line-segments  on  a  square  lattice.  We  have  also  t>bserved  the 
jamming  coverage  to  decrease  with  the  chain-length,  however,  we  concentrate 
here  primarily  on  percolation  quantities  for  monodisperse  chains.  In  the 
following  sectiim  we  describe  the  model.  In  section  three  the  results  of  the 
simulaticm  arc  presented  with  a  summarv  and  discussion  in  section  linn . 


2.  Model 

We  consider  a  square  lattice,  and  sticks  (i.c.  stiff-chains)  as  the  percolating 
units.  In  classical  percolation  theory  one  studies  the  statistics  of  a  single 
percolating  site.  We  wish  to  study  the  statistics  of  percolating  chains  u  here  a 
chain  is  defined  as  multiple  percolating  sites  connected  in  a  linear  fashion.  Thus 
the  length  ot  a  chain  is  the  number  of  linearly  connected  sites  on  the  lattice: 

=  1  corresponds  to  standard  site  percolation.  For  a  Hxed  chain-length  L  .  we 
generate  clusters  ot  connected  chains  in  the  following  way.  The  decision  to 
place  a  chain  is  determined  with  a  probability  p  with  the  help  of  pseudo¬ 
random  numbers.  It  the  landom  number  is  less  than  or  equal  to  p.  then  we 
randomly  select  one  of  the  tour  directional  orientations.  The  chain  is  then 
placed  along  the  chosen  direction  if  no  occupied  lattice  site  is  in  the  wav  as  two 
chains  are  not  allowed  to  share  a  lattice  site.  If  the  randomly  selected 
orientational  direction  fails,  then  we  choose  another  empty  site  randomly,  and 
a  random  orientational  direction.  Therefore,  a  chain  is  not  placed  if  it  were  to 
intersect  with  a  previously  occupied  lattice  site.  This  process  of  selecting  an 
t-iiipty  site  with  probability  p.  and  one  of  the  orientational  directions  randomly, 
and  attempt  to  place  the  chain  is  repeated  again  and  again,  for  a  lixed 
concentration  of  occupied  sites.  The  consequence  ot  this  rule  i:'.  that  a  jamming 
concentration  limit  p^  is  reached  when  it  is  not  possible  to  place  additional 
chains  in  the  lattice.  The  jamming  coverage  is  less  than  one  for  /.  larger  than 
one.  We  concentrate  here  mainly  on  the  percolation  quantities,  and  focus  on 
the  cluster  size,  their  number,  percolation  probability  and  the  second  moment 
of  the  cluster  size  distribution.  In  this  study,  a  cluster  is  defined  by  joining  the 
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nearest-neighbor  occupied  sites,  thus  two  neighboring  sticks  will  be  part  ot  the 
same  cluster  if  they  are  separated  by  one  lattice  constant. 


3.  Results 

Simulations  are  carried  t)ut  on  CONVEX  C1XP2  machine  with  32  MB  of 
memory.  This  caimputcr  ct>ntains  a  vector  processor,  yet  we  were  unable  to 
make  ma.ximum  use  due  to  the  hick  of  vector  code  in  the  algorithm.  For  each 
probability  p.  two  hundred  independent  runs  are  used  on  a  lOOO  ^  KHM)  lattice 
with  various  size  chtiins.  A  typical  computation  requires  two  hours  of  process¬ 
ing  time  for  each  set  of  realizations.  Twenty-eight  different  values  of  p  were 
used  for  seven  different  chain-lengths.  About  400  hours  of  computer  time  was 
used  in  this  project. 

As  we  mentioned  above,  we  have  analysed  percolation  probability  /’  . 
cluster  number,  and  the  second  moment  S={.s')  of  the  cluster  size  dis¬ 
tribution.  The  percolation  threshold  was  determined  by  hnding  the  peak  in  the 
second  moment.  The  result  of  the  percolation  threshold  versus  chain-length  is 
shown  in  Hg.  1.  We  see  that  the  percolation  threshold  decreases  on  increasing 
the  chain-length.  When  shown  on  a  log-log  plot  this  decay  exhibits  a  exponent 
power-law  dependence  i.e.  '.  This  exponent  .v  was  found  using  a 

least-squares  technique  to  be  about  one  half. 

The  exponent  beta  is  determined  by  examining  the  percolation  probability 
data  with  P,  —  {  p  ~/\  )^.  Fig.  2  shows  this  percolation  probability  i.e.  prob¬ 
ability  that  a  site  is  a  member  of  the  infinite  cluster  versus  the  probability  that  a 
site  is  occupied  for  various  chain-lengths.  Again,  it  can  be  seen  that  the 
percolation  threshold  decreases  with  increasing  chain-length.  It  should  also  be 
noted  that  the  percolation  probability  saturates  at  a  value  of  p  less  than  one. 
The  saturation  value  of  P .  is  the  jamming  coverage  jb],  and  it  decreases  with 
increasing  chain-length  as  expected. 

From  the  least-square  tit  of  the  log-log  plot  of  P,  versus  { p  -  p^.)  we 
estimate  (5.  We  have  evaluated  ft  for  various  chain-lengths  and  the  results  are 
presented  in  table  1.  Within  the  range  of  the  deviations,  these  values  for  fS 
seem  to  remain  the  same  for  different  chain-lengths,  and  is  consistent  with  the 
accepted  values  for  site  and  bond  percolation  (IJ. 

Variation  of  the  second  moment  S  (i.e.  the  susceptibility)  with  concentration 
p  is  shown  in  fig.  3  for  various  chain-lengths.  Tt^e  concentration  where  S 
reaches  its  maximum  value  is  the  percolation  threshold.  Again,  it  is  clear  from 
the  data  that  the  percolation  threshold  decreases  with  increasing  chain-length. 
A  log-log  plot  of  S  versus  {p  -  p^)  gives  the  exponent  y,  S  —  { p  -  p^)  The 
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Fig.  1.  Percolation  threshold  p^  versus  chain-length  (a)  linear  scale,  and  (b)  log-log  scale. 


values  of  y  for  various  chain-lengths  are  presented  in  table  1  which  appears  to 
show  no  change  with  chain  length  ( +  and  -  denote  the  estimates  at  p  above 
and  below  respectively). 

Using  the  scaling  relations  (I|,  other  exponents  like  the  correlation  length 
exponent  p  can  oe  estimated.  The  fractal  dimensionality  D  may  also  be 
evaluated  from  D  =  d  -  p/p.  The  computed  values  for  p  and  D  for  various 
chain-lengths  are  collected  in  table  1  for  various  chain-lengths. 
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4.  Summary 

We  have  presented  a  computer-simulation  study  of  the  percolation  ol  stiff 
chains/ sticks  in  a  two-dimensional  lattice.  The  percolation  exponents  seem  to 
remain  unchanged.  The  percolation  threshold,  on  the  other  hand,  depends 
strongly  on  the  chain-length  with  a  power-law  p^  —  '  .  The  jamming 

coverage  decreases  with  the  chain-length.  We  plan  to  study  the  percolation  of 
random  chains  in  the  near  future. 
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